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Abstract 
Introduction. The development of the polar regions of the World Ocean contributed to an increased interest in studying 


wave processes in water bodies with ice cover caused by the action of a mobile load. In most papers of domestic and 
foreign scientists, the ice sheet was considered as an elastic or viscoelastic plate loaded with a rectilinearly moving 
vertical force. However, when modeling the impact of vehicles on the ice cover, it is of interest to investigate problems 
in which the force moves along a more complex trajectory. Therefore, this study aims at developing a method for 
studying the behavior of the ice cover under the action of a force moving along a trajectory of a complex shape, obeying 
an arbitrary law of motion. 

Materials and Methods. A method for solving problems of the action of an arbitrarily moving force on the ice cover of 
a reservoir of finite depth filled with an inviscid incompressible fluid is proposed. The ice cover was considered as a 
viscoelastic plate lying on the surface of a liquid in a state of potential flow. A concentrated force moving along an 
arbitrary closed trajectory and being periodic in time was applied to the upper surface of the plate. Hydrodynamic 
pressure acted on the lower surface of the plate from the liquid side. Due to the periodicity of the load applied to the 
plate, an integral time transformation was used to solve differential equations describing the behavior of this system. 
Further, using traditional methods, formulas were obtained for calculating stresses and displacements in the plate and 
components of the velocity vector of liquid particles. These formulas were presented in the form of an iterated integral. 
Numerical methods were used to calculate integrals. 

Results. Calculations showed that the deflection of the ice cover increased markedly with the growth of speed and 
tangential acceleration of the load movement. An increase in the relaxation time of ice and a decrease in the radius of 
the trajectory of the load also caused an increase in deflection. The distribution of the fluid particle velocity vector over 
the depth of the reservoir was calculated. 

Discussion and Conclusions. The proposed method has shown its efficiency in solving problems about the impact of a 
moving load on the ice cover of a reservoir. With its help, the influence of the curvature of the trajectory of motion and 
the mechanical properties of ice, the kinematic characteristics of the movement of the load on the deflection of the ice 
cover was investigated. The method simulates the impact of vehicles on the ice cover of a reservoir. The results of its 


application can be used under the construction of ice roads or airdromes on ice. 
Keywords: infinite ice cover, moving load, arbitrary closed trajectory, variable speed. 
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AHHOTauHA 

Beedenue. Ocpoenve MonApHElx paltioHop Muposoro oKeaHa CMOcOOcTBOBAaIO HOBbIINeCHHIO MHTepeca K M3y4eHHIO 
BOJIHOBBIX IIPOIeCcCOB B BOJOeMAaX C JICXAHbIM MIOKPOBOM, OOYCJIOBJICHHBIX JelCTBHeM MOJBYWXKHOM Harpy3KH. B 
OOJIbIIMHCTBe paOOT OTeYECTBEHHBIX HW 3apyOexKHbIX YYCHBIX JIeqAHOM MOKPOB paccCMaTpHBaJIcA KaK ylipyrad WIM 
BA3KOYUpyrad WJiacTHHa, HarpyKeHHat MpAMONMHeMHO BWKyWelicd BepTHKaIbHOM cHoM. Ognako mp 
MOJCJIMpOBaHHH BO3eHCTBUA TpaHCIOpTHBIX CpeACTB Ha JIeqAHOM MOKpOB MpecTaBlaeT MHTepec paccmMoTpeHne 
3aa4, B KOTOPBIX CHJIa JBMKeTCA 10 Oosee CO%KHOM TpaekTopun. IlosToMy WebIO JaHHOrO MCCeAOBaHUA ABJIACTCA 
pa3pa0oTKa MeTOJa UccueqOBaHHA NOBeeHHA JeqAHOTO MOKpoOBa MOA WelMCTBHEM CHJIbI, JBMOKyIUelca 10 TpaeKTOpHu 
CIOXKHOM POPMBI, HOLYHHAACh MPOVM3BOJIbBHOMY 3aKOHY JBYDKeHHA. 

Mamepuanoit u memoooi. Upeqioxer MetToy pelieHua 3aa4 O TelCTBHM JBYKYLWIelicaA TIPOM3BOJIbHBIM OOpa30M CHIIBI 110 
JI€IAHOMY MOKpOBy BOJOeCMa KOHCYHOM TIyOMHBI, 3allOJIHeHHOTO HeBASKOM HecxKHMaeMOl 2K AKOCTHIO. JleqAHOM MOKpoB 
paccmMaTpHBalicd KaK BA3KOYIIpyrad@ IWiacTHHa, JlexKalllat Ha MOBEPXHOCTH 2 KHKOCTH, HaxOsAllelica B COCTOAHHM 
MOTeHIMaIbHoro TeyeHit. Ha BepXHIOIO NOBEPXHOCTb IWJIaCTHHbI MpHIO%KeHa COCpesOTOYeHHa CHa, TBWOKyWaica 10 
MIPOH3BOJIbHOM 3AMKHYTOM TpacKTOPHH VM ABIIAIOWMIAACA MepHOAMYeCcKON Ho BpemMeHu. Ha HYKHIOIO HMOBEPXHOCTb TIaCTHHbI 
CO CTOPOHbI %KHAKOCTH JelicTByeT THqpoqMHamuyecKoe JaBsieHHe. B cusly MepHOAM4YHOCTH NpHNOKeHHOM K TWIacTHHe 
Harpy3kKH TPH pelieHun udpdepeHuMabHbIX ypaBHeHHi, OMMCbIBaIOWIHX MOBeyeHHe WaHHOM CHCTeMbI, MCTIOJIb30BaJIN 
WHTerpasibHoe ImpeoOpa3o0BaHve m0 BpeMeHH. JJaee, NPHMeHAA TpawMUMOHHbIe MeTOJbI, NOMyYMIM COpMybl [yA 
BBIYHCIICHHA HallpsKeHHM HW MepeMeLeHHH B IWIaCTHHe H KOMIOHEHT BeEKTOpa CKOPOCTH YacTHL %KHAKOCTH. ITU PopMysBI 
TIpecTaBMIM B Be MOBTOpHOro MHTerpaya. Jit BEIYHCIeCHHA MHTerpasiOB MCHOMb3OBaJIN YMCICHHbIC METOJEL. 
Pezyismamot uccredosanua. Pacuetel MoKa3asIM, YTO MporHO WeqAHOrO MOKpOBa 3AMeTHO YBeIMYMBaeTCA C POCTOM 
CKOpOCTH UM KacaTeIbHOrO yCKOpeHHA JBHKeHHA Harpy3KH. K pocty nporuba TakoKe IPUBOAAT yBeMyeHHe BpeMeHH 
pewlakcallMM jIbya MW yMeHbIeHHe payuyca TpaeKTOpHu ABwKeHHA Harpy3KH. PaccuuTaHo paciipeyeseHve BeKTOpa 
CKOPOCTH 4acTHL *XHAKOCTH M0 rryOuHe BosOema. 

Oécystcoenue u 3akio4uenua. I[penioxKeHHbIM MeTOA MoKa3aI CBOIO 9*pPeKTMBHOCTh MPH pelwieHHM 3ayjay oO TelicTBHH 
NOABWKHOM Harpy3KH Ha JIeqaHOM MOKpoB BofOemMa. C ero MOMOLIbIO HCCIeqOBaHO BIMAHHe Ha UporHMO ezAHOrO 
MIOKPOBa KPHBH3HbI TpacKTOpHH ABHKeHHA M MeXaHHYeCKHX CBOMCTB JIba, KHHEMaTHYCCKHX XapakTepHCTHK JBWOKeHHA 
Harpy3KH. Metoy, MojeMpyeT Bo3yelcTBHe TpaHCHOpTHBIX CpeJICTB Ha JIeqAHOM MOKpoB BoyOema. Pe3ybTaTHI ero 


TIpHMeHeHHA MOryT ObITb HCHOJIb30BaHbI TIPH CTPOHTCJIBCTBE JICQOBBIX JOPOr WJM aSpPOAPOMOB Ha JIb]Ty. 


KoroueBbie c1I0Ba: OeCKOHeUHBIN Je AHOU TIOKPOB, ABYMKYUWaACA Harpy3Ka, MpOW3BOJIbHaA 3aMKHYTaA TpaecKTOpHA, 


TlepeMeHHawA CKOPOCTb. 


BuaarogapHocrn. ABTOp BbIpakKaeT OnarojyjapHOcTb peueH3eHTaM 3a YKa3aHHble 3aMCYaHHA, KOTOPbIe MO3BOJIMJIM 


TLIOBbICHTb KAYCCTBO CTaTbH. 


Aaa unTupopanua. Tanadypoun A.B. Metoy peulenua 3aqauu 0 JBWKeHHM Harpy3KH MO jeqtHOMy MOKpoBy 
BoqoOeMa 0 CiIOKHOM Tpaextopuu. Advanced Engineering Research (Rostov-on-Don). 2023;23(1):34-40. 
https://doi.org/10.23947/2687-1653-2023-23-1-34-40 


Introduction. Recently, in connection with the development of the northern territories and areas of the World 
Ocean, great attention has been paid to the studying the behavior of the ice cover under the action of external load. A 
large number of works by domestic and foreign scientists are devoted to these issues. The impact of a pulsed moving 
load on a viscoelastic floating plate was considered in [1]. The effect of a moving load on the ice cover frozen to the 


channel walls was presented in [2, 3]. The fluctuations of the ice cover caused by a load moving at a constant speed 
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were considered in [4]. The behavior of a semi-infinite ice sheet under the action of a uniformly moving load on it was 
studied in [5]. The propagation of waves excited along a channel with an ice cover was discussed in [6]. In [7], 
nonlinear models were used in simulation. The monograph [8] presented the results of studies of surface waves in a sea 
with a floating broken and solid ice cover. The ice failure under the action of a moving load was considered in [9]. In 
the above works, the rectilinear load movement was mainly considered. The objective of this study was to develop a 
method for solving the problem of the impact of a force moving across the ice cover in an arbitrary way. Achieving this 
goal makes it possible to more accurately simulate the impact on the ice cover of vehicles, whose movement often 
occurs along rather complex trajectories and according to a complex law. 

Problem Statement. The oscillations of an infinite ice cover lying on the surface of a reservoir of finite depth under 
the action of a force moving arbitrarily along a closed trajectory were considered. The ice cover was modeled by a thin 
viscoelastic plate, whose mechanical properties were described by the Kelvin-Voigt model. The reservoir was filled 
with an incompressible liquid. 

Materials and Methods. The ice cover bending was described by the differential equation [7]: 

(1+ ,0;)A?W + c?ZW + kW + bd,®|,-9 =”, (1) 
where W(x,y,t) — ice cover deflection; D=Eh'{ 12( 1-1”) ); E — Young's modulus; 1 — Poisson's ratio; h — thickness 
of the ice cover; t, — relaxation time of deformations; A?=0,’+20,7 a7 +0,*; P,— ice density; p, — water density; c= 
p,h/D, k=p,.g/D, b= p,/D; P(x,y,t) — load acting on the ice surface; ®(x,y,z,t) — potential of fluid movement. The fluid 
behavior was subject to the equation: 

A®=0. 

Boundary conditions at the ice-water interface at z=0 and at the bottom of the reservoir z=—H (H — reservoir depth) 

had the form: 
0,W=0,Q|.-o, 0-D| ._., =0. 

It was supposed that force P (x, y, t) moved along an arbitrary closed trajectory y in an arbitrary way. It was assumed 
that P=P(s(t)), where s — arc coordinate measured from some fixed point of the curve y. If function s(t) was periodic 
with period T, then it was quite obvious that value P(s(t)), would also be periodic, and also with period T. The 
parametric assignment of the trajectory was taken as: 

x = X(t) 
t = Yo(t) 
where ¢ — time. 

Considering the steady-state process and applying the integral transformation with respect to variables x and y as 
well as the final integral transformation with respect to t on the interval [0; 7], we obtained: 

(1— i@)t) p'Wo— °°,” Wot kKWotbiaj,Do|--c=Po/D, 
8/@Mp — p’ Do=0, 
where Wy, ®) — images of unknown functions W and ®, p*=’+ a’; 4; a — parameters of integral transformations, 
respectively, by variables x and y; o,=2an/T, n=0, 1, 2,3, ... . Solving the second equation allowing for the boundary 
conditions, we obtained: 
Do=-i0,,Woch(p(z+H))/psh(pH). 

Then, the first equation followed: 

Wo=PA((c°+beth(pH Vp) ©, -p*-k+i OnTop’). 
The moving concentrated force was approximated by the function: 

P(x,y,t)=€° exp(-€°((x-X,(t)) + (y-yolt))”))/a 
where ¢ — parameter. 

Further, by performing quite obvious transformations, we got: 


1 i 2 ( 7,,cos(@,(t—-T)) %,Sin(@, (t—T seid 
w= sofll esl ( ( ve 2 ( ( 2, (ome Jos 


n=l r, 


n n 


rin=(c+bcth (PHYp) Oy P'-k,, ry= top", R'=5 +B, 
5=x(T-x, B= yo(t)V. 
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Knowing W, it was possible to determine formulas for calculating the components of the displacement vector and 
the stress tensor at any point of the ice cover using available ratios. 

In the numerical implementation of the proposed method, a problem arose related to the need to calculate integrals. 
One of them as a subintegral function had a strongly oscillating function for sufficiently large n, and the second was 
improper with an infinite upper limit. 

To calculate integrals from strongly oscillating functions, a quadrature formula was used, obtained by the cubic 


spline method [11]: 


b ; 1 N-1 eon =p 

10x 
Jef (x)dx~- - (M,.-M,). 
a @ j=l Si 


where h; — lengths of elementary segments into which the interval was divided [a; b]; S(x) — approximation of f(x) by 
cubic spline, M; = 5 (Xp): 
When calculating the improper integral, approximate ratio ss f(p)dp = f f(p)dp was used for a sufficiently 


large A. Value A was chosen so that the estimate of the allowable error | if i. f (p)dp| was small enough. 

Using traditional methods, it was possible to obtain the required estimates for the calculated quantities. For example, 
for deflection value of the ice cover, this estimate had the form: 
(c? +bcth(AH)/A)o, +A*(1+0,1, )+k 


no 


loa) 


fw(p) dp 


A 


_A2s4p2 
<4e2e 4 


[(? +beth(AH)/ A)@, — A’ -k] +[o,1,A*} 


w=) w(p)dp + J,” w@)dp. 
When calculating the sum of a series, the Lanczos sigma multiplier method was used to accelerate its convergence. 
Research Results. The calculations were carried out for the case of the action of a single concentrated force that 
moved along a closed curve, shown in Figure 1: 


“si =f; 
cae 
2sin( Tp) 
It was assumed that the thickness of the ice cover was h=0.25 m, Young's modulus of the plate material 
E=500,000,000 N/m’, Poisson's ratio u=1/3, ice density p=900kg/m*, water density p=1,000 kg/m’, reservoir 


depth H=5 m, e=2.5. The radii determining the shape of the force trajectory were assumed to be R;=15 m, Ro=9 m, 
R3=3 m (Fig. 1). 


ans, te[0:T]. 


Y,m 


4 -3 -2 -1 O 1 2 3 4 X,m 
Fig. 1. Trajectory of the concentrated force 


Figure 2 shows the change in the deflection of the ice cover (a), as well as the maximum values of normal stresses 


(at z=th) Sx and Sy (b and c) and at the movement speed of single concentrated force v =7.4022 m/s and tangential 


acceleration w,=0. 


Mechanics 


http://vestnik-donstu.ru 


Advanced Engineering Research (Rostov-on-Don). 2023 ;23(1):34—40. eISSN 2687-1653 


Sx, N/m? Sy, N/m’ 
, } | 
200 | ’ 
200 \ o| | 
iy -200 1 i iy 
” -200 “6 
ioe Y,m ilies 
-50 aie ; — O* 
ae E se () Ne Es ' 
~~ = ei Y,m 
Xin -50 X,m~-50 0 
a) b) c) 


Fig. 2. Change in deflection and stresses of the ice cover: a — change in deflection; 
b — change in stress S,; c — change in stress S, 


On these graphs, the red dot indicates the place of application of force, and the blue color indicates the trajectory of 
movement. 

Figure 3 shows the dependence of the change in the maximum value of the deflection of the ice cover on the force 
speed rate (a) at the moment t=T/2. The position of the force on the trajectory at the moment is marked by a red 
dot (Fig. 1), while the tangential acceleration w,=0. Two graphs are given, one of which corresponds to the radius of the 
trajectory R,=15 m (solid line), and the second — to radius R;=5 m (dotted line). Graph b in Figure 3 shows the 
dependence of the maximum value of the ice cover deflection on the tangential acceleration rate at moment t=T at the 
point of the trajectory (2; 0) in Figure | at speed v=0. The graph in Figure 3 shows the dependence of the maximum 


deflection value on the value of the strain relaxation time. 
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Fig. 3. Dependences of the ice cover deflection on: a — speed; 
b — tangential acceleration; c — relaxation time 


The liquid behavior is shown in Figure 4. It illustrates the distribution of the speed vector of liquid particles over the 


depth of the reservoir at t=T/2, force speed v=3.7011 m/s and acceleration w,=0. 
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Fig. 4. Distribution of the speed vector of liquid particles over the reservoir depth 


Discussion and Conclusions. The influence of the movement trajectory curvature and the mechanical properties of 
ice on the ice cover deflection, speed and acceleration of the load movement was investigated. Calculations have shown 
that the deflection of the ice cover depends significantly on the movement speed of the force and acceleration of 
movement. 

However, the nature of qualitative changes in displacements and stresses in the ice cover caused by the action of a 
moving force with a change in speed and acceleration of movement changed slightly. 

Mechanical properties of ice, specifically, the relaxation time, had a noticeable effect on the deflection of the ice 
cover. 

The results obtained and the proposed method for solving such problems can be used in the construction of ice roads 
or airfields on ice. Moreover, the proposed solution method has shown its effectiveness and can be used to solve other 


similar problems. 
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